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ABSTRACT 

Deep observations of the Small Magellanic Cloud (SMC) and region were carried 
out in the hard X-ray band by the INTEGRAL observatory in 2008-2009. The field 
of view of the instrument permitted simultaneous coverage of the entire SMC and the 
eastern end of the Magellanic Bridge. In total, INTEGRAL detected seven sources in 
the SMC and five in the Magellanic Bridge; the majority of the sources were previously 
unknown systems. Several of the new sources were detected undergoing bright X- 
ray outbursts and all the sources exhibited transient behaviour except the supergiant 
system SMC X-l. They are all thought to be High Mass X-ray Binary (HMXB) systems 
in which the compact object is a neutron star. 
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1 INTRODUCTION AND BACKGROUND 

The Small Magellanic Cloud (SMC) present a nearby, easily 
observable irregular galaxy in which it is possible to study 
the results of tidal interactions of star formation. Multiwave- 
length studies of the SMC have shown that it contains a 
large number of X-ray binary pulsars all but one of them 
in Be/X-ray binary systems (Coe et al. 2005, 2009). This 
preponderance of such systems is most likely explained by 
tidally-triggered star formation induced by the most recent 
close approach of the SMC and LMC (Gardiner & Noguchi 
1996). However, the most recent close approach of the SMC 
and LMC was 200 Myr ago - much longer than the evo- 
lutionary timescale of Be/X-ray binaries. This implies that 
either there has been a significant delay between the en- 
counter of the SMC and LMC and the onset of star forma- 
tion, or that subsequent waves of star formation have given 
rise to these Be/X-ray binaries (Harris & Zaritsky 2004). 



These objects, being tracers of star formation (Grimm et al. 
2003) , give direct insights into the star formation history of 
their host galaxies. 

In this paper an extended set of deep observations of 
the SMC region are reported. One major new discovery from 
this INTEGRAL study was a population of HMXBs in the 
Magellanic Bridge. These observations have been published 
elsewhere(McBride et al, 2010). This paper focusses exclu- 
sively on the HMXB population of the SMC itself. 



2 OBSERVATIONS 

The IBIS telescope (Ubertini et al. 2003) on board INTE- 
GRAL, which is optimised for an energy range of 15-200 
keV and has a field of view of 29x29 degrees (Full Width 
Zero Intensity response), is well suited to observing large 
sky areas for point sources. As part of a key programme 
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Table 1. INTEGRAL observation log. 



Rev Number 


MJD start 


Exposure (ks) 


745 


54788.7 


43.8 


746 


54791.8 


47.6 


747 


i y±. o 


47.2 


748 


54797.7 


46.8 


749 


J40UU.O 


48.0 


750 


^480^ 7 


55.7 


751 


^4S0fi 7 

JIOUU . t 


64.6 


752 




57.7 


753 


54812.7 


58.8 


751 


54815.7 


45.3 


755 


54818.7 


60.9 


756 


54821.7 


57.2 


796 


54941.4 


79.6 


797 


54944.4 


62.7 


812 


54989.2 


63.3 


813 


54992.2 


53.1 


814 


54995.2 


60.8 


815 


54998.2 


58.8 


816 


55001.2 


62.1 


817 


55004.2 


63.1 


818 


55007.2 


40.2 



Table 2. Table of definite INTEGRAL source detections in the 
SMC 



Source name Times of detection (MJD) 



SMC X-l 


Various - see Figure [T] 


SXP11.5 


54989 - 55007 


SXP756 


54941 - 54947 


IGR J00515-7328 


54800 - 55000 


SXP6.85 


54790 - 54820 



monitoring campaign on the SMC and 47 Tuc, INTEGRAL 
observed the SMC and Magellanic Bridge for approximately 
90 ks per satellite revolution (~3 days) from 2008 November 
11 to 2009 June 25 - see Table [T] for a detailed journal of all 
the observation dates and exposures. 

Individual pointings (science windows) were processed 
using the INTEGRAL Offline Science Analysis v.7.0 (OSA, 
Goldwurm et al. 2003) and were mosaiced using the weighted 
mean of the flux in the 3-10 keV (JEM-X) and 15-35 keV 
(IBIS) energy ranges. Proprietary software was used to mo- 
saic the observations from successive revolutions to improve 
the exposure and thereby the sensitivity to faint sources. 
Lightcurves in these energy bands were generated on sci- 
ence window ( 2000 s) and revolution time-scales. The IBIS 
energy band was chosen to maximise the detection signifi- 
cance of SMC X-l, and hence other SMC accreting X-ray 
pulsars, which have similar spectral shapes to SMC X-l in 
this energy range. An IBIS mosaic of data from revolutions 
752-756 in the 15-35 keV is shown in McBride et al (2010). 



3 INDIVIDUAL SOURCES SEEN IN THE SMC 

A list of all the definite source detections and when they 
occurred is presented in Table [2] 



3.1 SMC X-l 

As a neutron star accreting from a supergiant overfilling its 
Roche lobe, SMC X-l is the only persistent accreting X- 
ray pulsar in the SMC. This transfer of angular momentum 
serves to spin up the pulsar, and SMC X-l has demonstrated 
consistent spin up over the last 40 years. It has a 3.89 d 
orbital period with eclipses lasting 1.12 d (Wojdowski et al. 
1998) and also shows superorbital modulation at ~ 60 d, 
thought to be the result of a precessing, warped accretion 
disk (Clarkson et al. 2003). 

The first set of IBIS observations cover almost a full 
superorbital cycle, while the subsequent observations occur 
predominantly in the superorbital low state. See Fig. Q] for 
the whole lightcurve and Figure [2] for details of the binary 
eclipse. 

Folding data from revolutions 745 to 749 on the orbital 
ephemeris from Wojdowski et al. (1998) resulted in an or- 
bital profile during superorbital high state. (See Figure [2] 
top panel). For comparison, an orbital profile during super- 
orbital low state was generated using data from revolutions 
751 to 756 and 812 to 818. There are significant differences 
between the eclipse profiles during these two superorbital 
phases. Eclipse ingress and egress are steep, almost square, 
during superorbital high state, while they are much more 
gradual (almost 0.1 phase longer) during superorbital low 
state. These results are in agreement with the orbital profile 
as a function of superorbital phase presented in Trowbridge, 
Nowak & Wilms (2007). These authors use softer X-ray data 
(0.2-12 keV) to show that the eclipse is wide and shallow dur- 
ing superorbital low state while it is narrow and deep during 
superorbital high state. 

From these data there is no evidence that the eclipse in 
hard X-rays (15-35 keV) is shorter than in soft X-rays (0.2- 
12keV, Trowbridge, Nowak & Wilms 2007), which may be 
expected due to the increased transparency of the supergiant 
photosphere to hard X-ray photons. There may, however, be 
some evidence for asymmetry of ingress and egress during 
the high superorbital state eclipse. The eclipse profile in Fig- 
ure [2] shows a similar shape to 20-40 keV eclipse profile in 
McBride et al. (2007). 

An IBIS observation at MJD 54797, during the super- 
orbital bright phase, was used to determine the current spin 
period of the pulsar in SMC X-l. This was determined using 
an epoch folding technique to be 0.702010 ± 0.000008 s, and 
is plotted above the arrow in Fig [3] A straight line model 
does not fit the period evolution very well, but gives a good 
idea of the general spin up of the system over its lifetime, 
which can be estimated as ~ —3.8 x 10 -4 s/yr. 



3.2 SXP756 

A source was detected in INTEGRAL at a position coin- 
cident with SXP756 (Coe & Edge, 2004) within an uncer- 
tainty of ~1 arcminute over two satellite orbits (MJD54941 - 
54947). Follow-up RXTE observations (MJD54943) revealed 
evidence for weak pulsations around 730±5s, the relatively 
large uncertainty reflecting the fact that the observation ex- 
tended over 13ks and therefore does not cover many period 
cycles. Figure [3] shows the Lomb-Scargle periodogram for 
the RXTE/PCA data. 

Independent confirmation that both the INTEGRAL 
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Figure 1. The IBIS measurements of SMC X-l (diamonds). 
Each measurement represents the flux (15-35 keV) averaged over 
a spacecraft orbit. The RXTE/ASM lightcurve, arbitrarily scaled 
is overplotted in grey to show the superorbital modulation and 
the binary eclipses. 



Figure 3. Upper panel: The history of the spin period of SMC 
X-l. The latest IBIS measurement on MJD 54797 is indicated 
by an arrow. The penultimate measurement is from the previous 
INTEGRAL observations in 2003 (McBride et al, 2007) Lower 
panel: the deviations of the period from a simple linear trend. 
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Figure 2. . IBIS 15-35 keV data of SMC X-l folded on the orbital 
period during superorbital high state (top panel) and superorbital 
low state (lower panel). The binning is 0.05 of the orbital period. 



and the RXTE sources are probably SXP756 comes from 
the OGLE III optical monitoring of the counterpart. This 
source is unusual in that it shows striking optical behaviour 
around the time of the periastron passage of the neutron star 
(Coe & Edge, 2004). These authors quote an ephemeris for 
outbursts of T = JD2449830 + 394N. The outburst reported 
here peaked at JD2454945 which gives a phase for the INTE- 
GRAL observations of 0.98 using the published ephemeris; 
coinciding with the sharp optical flare seen in the OGLE III 
data at that time - see Figure [5] Thus there can be no doubt 
that it is SXP756 that was seen in the INTEGRAL data. 

The significant neutron star spin up in this system is 
shown from our long-term RXTE monitoring programme - 
see Figure [6] Measurements up to MJD54000 come from 
Galache et al 2008. From the data shown in this figure it is 
possible to determine a straight line fit using a x 2 routine 
and establish an average P = -1.28s/year over ~11 years. 




0.002 0.004 0.006 0.008 

Frequency (1 /s) 

Figure 4. Lomb-Scargle periodogram of RXTE PCA data for 
SXP756 covering the period range 100-1000s. The highest peak 
corresponds to a period of 730±5s. 



Using Figure 5 of Coe, McBride & Corbet (2010) this value 
implies a typical outburst P of ~5 times that value. 

The X-ray luminosity was calculated from the observed 
peak pulse amplitude in counts/pcu/s that occurred during 
the current outburst. This amplitude is converted to lumi- 
nosity assuming that 1 RXTE count/pcu/s = 0.4xl0 37 erg/s 
at the distance of 60kpc to the SMC (though the depth of 
the sources within the SMC is unknown and may affect this 
distance by up to ±10kpc - Laney & Stobie, 1986). The X- 
ray spectrum was assumed to be a power law with a photon 
index = 1.5 and an Nh = 6 x W 20 cm~ 2 . Furthermore it was 
assumed that there was an average pulse fraction of 33% for 
the system and hence the correct total flux is 3 times the 
pulsed component. From this we find that the approximate 
luminosity from our RXTE data in the 3-10keV range is 
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Figure 5. OGLE III I-band lightcurve of the optical counterpart 
to SXP756. Optical flares are clearly visible every 394d (Coe &: 
Edge, 2004). The INTEGRAL detection was coincident with the 
last optical flare seen in these data. 
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Figure 6. The pulse period history of SXP756. All the mea- 
surements shown, except the two marked with a triangle, are 
from RXTE. The triangle data points are from XMM (Habcrl 
& Pietsch 2004, Haberl, Eger & Pietsch 2008). The final point 
coincides with the INTEGRAL detection reported in this work. 



~ 2 x 10 37 erg/s. Using the harder IBIS data (15-35 keV) a 
similar luminosity of (~ 1) x 10 37 erg/s is found. 

Locating these values for P and L x during the outburst 
on the plot of logP versus logPL 3 ^ 7 for SMC sources (Figure 
6 of Coe, McBride & Corbet, 2010) puts SXP756 close to 
the end of the distribution in the region occupied by two 
other long period pulsars, SXP892 and SXP1323. 
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Figure 7. INTEGRAL/IBIS 16-35 keV light curve of IGR 
J00515-7328 



3.3 IGR J00515 7328 

A new source, IGR J00515-7328, was active when the third 
set of SMC observations started in INTEGRAL revolution 
812 (MJD54989). The source continued to be marginally de- 
tectable for 3 subsequent revolutions, peaking in flux around 
MJD54994 (June 12th 2009), see Figure [7] 

The combined significance across these 4 revolutions is 
about 5.9 sigma, a robust detection of a faint source. The 
IBIS source position is RA 00:51:28, Dec -73:30:33 with an 
error circle of radius 4.2 arcmin. The source had faded below 
the IBIS detection limit by MJD55000. 

A short 500s Swift/XRT TOO was performed on 
MJD55003, shortly after the source had faded below the de- 
tection threshold in IBIS. One faint source was seen within 
the IBIS error circle at RA=00:51:59, Dec=-73:29:25 with a 
5.2" radius error circle. With detection based on 40 pho- 
tons, further diagnostic spectral and/or timing analysis was 
infeasible. A basic extrapolation of the IBIS and XRT fluxes, 
assuming a typical SMC Be-X spectrum (power law with a 
photon index of -1.0) and a steady decline of the flux from 
the source after the IBIS detection, shows the two fluxes to 
be compatible, but we cannot unambiguously associate the 
IBIS and XRT detections. The XRT source is very likely as- 
sociated with a ROSAT PSPC source (Haberl et al., 2000) 
and a ROSAT HRI source at RA=00:51:59.6 Dec=-73:29:29 
(Sasaki et al., 2000), see Figure [8] The flux measured by 



HRI (~ 5 x 10 erg/s) is around 20 times lower than the 
XRT flux (10 36 erg/s), implying that while the XRT de- 
tection was during a decaying outburst, the HRI detection 
probably represents the quiescent flux from the source. 

It is also worth noting that there are 6 XMM sources 
in the IBIS error circle taken from the XMM Serendipitous 
Source Catalog (2XMMi version) (Watson et al, 2009). How- 
ever, none lie close to the area shown in Figure [H] Further- 
more these XMM sources are restricted to a band across 
the IBIS error circle suggesting incomplete coverage of this 
region by XMM. 

2MASS J00520059-7329255 lies 5.6" away from the 
XRT position, marginally outside the 90% error circle and 
consistent with the ROSAT HRI & PSPC positions. This 
source is also known as [M2002] SMC 21983 (Massey 2002) 
and has B and V magnitudes of 15.18, typical of Be systems 
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Figure 8. A Digitised Sky Survey red-band image showing all 
the ROSAT and Swift sources known within the IBIS error circle 
for IGR J00515-7328. The position of 2MASS J00520059- 7329255 
is indicated by a white cross. The Swift/XRT error circle has a 
radius of 5.2" . 



period in the MJD54994 data. Again, there is no robust 
way to associate this periodicity with IGR J00515-7328. 

Overall, the detection of an X-ray source showing a high 
level of variability in flux between two epochs, and associ- 
ation with an optical counterpart which is typical of a Be 
star, suggests that IGR J00515-7328 is another high-mass 
X-ray binary system within the SMC. 



3.4 SXP6.85 

SXP6.85 (= XTE J0103-728) underwent a large Type II 
outburst beginning on 2008 August 10 and the source was 
consistently seen in RXTE observations for the following 20 
weeks (MJD54688 - 54830) and by INTEGRAL in the period 
MJD54790 - 54820. Townsend et al (2010a) presented X-ray 
timing and spectroscopic analysis from the Rossi X-ray Tim- 
ing Explorer (RXTE) and the INTEGRAL observatory. A 
comparison with the Optical Gravitational Lensing Experi- 
ment (OGLE) III light curve of the Be counterpart showed 
the X-ray outbursts from this source to coincide with times 
of optical maximum. 




2000 3000 4000 5000 

Date (MJD = JD 2450000) 

Figure 9. OGLE III I-band lightcurve of the proposed optical 
counterpart IGR J00515-7328 over ~8 years. 



in the SMC. The source to the NW of 2MASS J00520059- 
7329255 is also in the Mssey catalogue (no: 21891) and has 
V=16.47 & B-V=0.42 and is, again, a possible counterpart 
to a Be X-ray binary system. However, its location well 
clear of all three error circles makes its association with the 
ROSAT and Swift objects much less likely. The OGLE light 
curve (see Figure © of 2MASS J00520059- 7329255 (OGLE 
III name SMC 103.4 33693) shows considerable variability 
on several timescales. Despite this, timing analysis of the 
de-trended light curve yielded no significant periodicities. 

A search of XTE monitoring observations performed on 
MJD54988 (just before the IBIS detection) and MJD54994 
(during the outburst) found only a hint of a 8.29s pulse 



3.5 SXP11.5 = IGR J01054-7253 

A large outburst from a previously unknown source was 
seen by INTEGRAL starting MJD54989, lasting ~20d and 
with a peak 3-10 keV luminosity of L x = 7 x 10 37 erg/ s. A 
rapid follow-up observation using the Swift observatory re- 
fined the X-ray position and permitted the identification of 
an associated optical counterpart. Townsend et al (2010b) 
present X-ray and optical data on this proposed new Be/X- 
ray binary pulsar, IGR J01054-7253 = SXP11.5. The opti- 
cal counterpart was shown to be [M2002] SMC 59977 and 
OGLE III data on that source presented. Following the dis- 
covery, RXTE began monitoring SXP11.5 approximately 3 
time per week for the duration of the outburst (MJD = 
55001 - 55063). These data revealed a clear pulse period of 
11.48s (Corbet et al., 2009). In addition, a strong Doppler 
modulation of the 11.48s pulse period was seen and inter- 
preted as an accretion driven spin-up. Model fitting to these 
data allowed Townsend et al (2010b) to distinguish the spin- 
up rate from the orbital modulation and identify an orbital 
period of 36.3T0.4 days and an eccentricity of 0.28T0.03 
(Townsend et al., 2009). 



3.6 Tentative new sources 

IBIS 15-35 keV and JEM-X 3-10 keV maps on time-scales 
varying between single revolutions up to five consecutive 
revolutions were searched for excesses which could indicate 
potential new sources in the SMC. A list of candidates is sup- 
plied in Table [3] and mapped on to the SMC in Figure 1101 
With the large error circles on these objects, it is not yet 
possible to constrain the nature of these sources. However, 
their transient nature makes it likely that they are binary 
systems in the SMC. 

One of these candidate sources (IGR J00523— 7217) was 
detected by combining data over a ~ 10 day period in June 
2009 and is spatially consistent with SXP327 & SXP4.78, 
known Be/X-ray binaries in the SMC. The region of the 
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Table 3. Candidate sources in the SMC. Coordinates are in J2000.0. For IBIS detections the fluxes are given in the 15-35 keV range, 
while JEM-X fluxes are quoted in the 3-10 keV range. 



Name 




RA 


Dec 


Significance 


Error radius 


IBIS 


JEM-X 


Flux 






hms 


dms 




arcmin (90%) 


MJD 


MJD 


(cts/s) 


ICR J00523- 


7217 


00 52 19 


-72 17 06 


5.3tr 


4.7 


54989-55001 




0.39±0.07 


IGR J01134- 


-7325 


01 13 06 


-73 25 17 


3.3ct 


1.8 




54812-54822 


(1±0.3) x 10~ 2 
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Figure 10. Sources detected by INTEGRAL overplotted on an 
SMC H I column density map (Stanimirovic, Staveley-Smith & 
Jones 2004). The sources in red are definite source detections 
discussed in this paper. The two sources in black are the candidate 
sources listed in Table [3] 



SMC in which this potential new source was seen was mon- 
itored weekly by RXTE, but pulsations were not detected 
from either of those sources. 



4 DISCUSSION 

The SMC is turning out to be an exciting nest of X-ray 
binary pulsars. It is possible to estimate the number of sys- 
tems one would expect based upon the relative masses of 
our galaxy and the SMC. This mass ratio is approx 50, so 
with 64 known or suspected systems in our galaxy we would 
only expect 1 or 2 systems in the SMC. However, Maeder, 
Grebel & Mermilliod (1999) have shown that the fraction of 
Be stars to B stars is 0.39 in the SMC compared with 0.16 in 
our galaxy. So this raises the expected number of Be/X-ray 
systems to approx 3 for the SMC - but we now know of more 
than 50 such systems in SMC! 

Many of these detections have come from XTE, Chan- 
dra & XMM X-ray observations over the last couple of years 
(Galache et al. 2008, McGowan et al. 2008, Haberl, Eger & 
Pietsch 2008). This large number suggests a dramatic phase 
of star birth in the past, probably associated with the most 
recent closest approach between the SMC and the LMC 
some 0.2 Gyrs ago (Gardiner & Noguchi 1996). Even more 
extreme, the very recent work of Naze et al (2003) in just 
one 20 x 20 arcmin Chandra field identified more than 20 
probable Be/X-ray binary systems. Multiplying these num- 



bers up by the 2x2 degree size of the SMC, and allowing for 
~10% X-ray duty cycles, suggests the final number of Be/X- 
ray binaries could be well in excess of 1,000! These observa- 
tions of the SMC are not only providing a great sample of 
HMXBs for study, but are also providing direct insights into 
the history of our neighbouring galaxy. Since the positional 
accuracy of Chandra & XMM is a few arcsec, or less, optical 
counterpart have been located for many of these potential 
XRB systems. 

The reason for the large number of HMXBs in the SMC 
probably lies in the history of the Magellanic Clouds. De- 
tailed HI mapping by Stavely-Smith et al (1997) and Put- 
man et al (1998) has shown a strong bridge of material be- 
tween the Magellanic Clouds and between them and our own 
galaxy. Furthermore, Stavely-Smith et al have demonstrated 
the existence of a large number of supernova remnants of a 
similar age (~5 Myr), strongly suggesting enhanced star- 
birth has taken place as a result of tidal interactions be- 
tween these component systems. It seems very likely that 
the previous closest approach of the SMC to the LMC ~100 
Myrs ago may have triggered the birth of many new mas- 
sive stars which have given rise to the current population of 
HMXBs. In fact, other authors (eg Popov et al 1998) claim 
that the presence of large numbers of HMXBs may be the 
best indication of starburst activity in a system. 

If the explanation of tidal interactions producing lots of 
starbirth is correct, it is very likely that there are many more 
systems waiting to be discovered in the SMC than the 60- 
100 known at this time. In total, INTEGRAL detected seven 
sources in the SMC - four of which were previously unknown 
systems. Combining this with the population of the five new 
sources detected in the Magellanic Bridge makes this a very 
productive survey of the SMC region. 

Comments on specific SMC sources are: 

• SMC X-l was overwhelmingly the strongest and most 
persistent source seen in this survey. It remains the only 
known supergiant system in the SMC and the results pre- 
sented here on the changing orbital profile strongly support 
the earlier work of Trowbridge, Nowak & Wilms (2007) at 
lower X-ray energies. It is clear that the superorbital pe- 
riod plays a major role in the obscuration of the emission 
from this system and the higher energy results presented 
here should provide further useful constraints on modeling 
this phenomenon. 

• SXP6.85 is typical of the transient sources detected in 
this survey. Townsend et al (2010a) attribute this flux vari- 
ability to the circumstellar disk increasing in size, causing 
mass accretion onto the neutron star. Ground based IR pho- 
tometry and Ha spectroscopy obtained during the outburst 
were used as a measure of the size of the circumstellar disk 
and lent support to this picture. In addition, the folded X- 
ray light curves seem to indicate complex changes in the ge- 
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ometry of the accretion regions on the surface of the neutron 
star, indicative of an inhomogeneous density distribution in 
the circumstellar material causing a variable accretion rate 
onto the neutron star. 

• The behaviour of SXP11.5 was similar to that of 
SXP6.85, with the addition of a clear Doppler period shift 
- something that is often obscured by much larger accretion 
torques during outbursts. The orbital and spin periods place 
this system nicely onto the Be/X-ray binary region of the 
Corbet diagram of binary pulsars. P is also consistent with 
that previously observed in other Be/X-ray binary pulsars 
in the Magellanic Clouds. 

The large number of sources seen by INTEGRAL is 
strongly indicative of the substantial advantage such a wide 
field instrument has when it comes to determining the ex- 
tent of a population of transients systems. Though parts of 
the SMC have been regularly monitored by RXTE over a 
decade (Galache et al., 2008), the smaller field of view, com- 
bined with the collimator response makes it difficult exactly 
to quantify the duty cycles of these HMXB systems. Other 
wide field instruments such as the Swift/BAT simply to do 
not have the combined spatial and flux sensitivity of INTE- 
GRAL/IBIS. Thus the conclusion of the INTEGRAL results 
presented here is that the total population of HMXBs in the 
SMC is probably much greater than the 60-70 systems so far 
catalogued, with consequent implications for the star forma- 
tion rate in this galaxy. Plus, there is probably a significant 
population along the Bridge - that may decline in number 
(along with the age of the stellar population) as one moves 
from the SMC to the LMC. Future INTEGRAL observa- 
tions should continue to enhance our understanding of all 
parts of the Magellanic system. 
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